Brans-Dicke Cosmology with a scalar field potential 
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Three solutions of the Brans-Dicke theory with a self-interacting quartic potential and perfect 
fluid distribution are presented for a spatially flat FRW geometry. The physical behavior is consistent 
with the recent cosmological scenario favored by type la supernova observations, indicating an 
accelerated expansion of the Universe. 
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Recent observations of type la supernovae with redshift up to about z < 1 provided evidence that we may live in 
"^j". a low mass-density Universe, with the contribution of the non-relativistic matter (baryonic plus dark) to the total 
energy density of the Universe of the order of f2 TO ~ 0.3 Q. The value of f2 m is significantly less than unity and, 
consequently, either the Universe is open or there is some additional energy density p sufficient to reach the value 
Qtotai = 1, predicted by inflationary theory. Observations also show that the deceleration parameter of the Universe 
q is in the range —1 < q < 0, and the present-day Universe undergoes an accelerated expansionary evolution. 

Several physical models have been proposed to give a consistent physical interpretation to these observational facts. 
Photons propagating in extra-galactic magnetic fields can oscillate into very light axions, so the supernovae appear 
dimer and more distant than they really are @. For the missing energy one candidate is the vacuum energy density, 
or the cosmological constant A [Q. A charged Universe, in which the long range force carriers have a small mass, 
has been considered as a cause for the cosmic acceleration in j5j. Another possibilities are cosmologies based on a 
mixture of cold dark matter and quintessence, a slowly- varying, spatially inhomogeneous component H. An example 
of implementation of the idea of quintessence is the suggestion that it is the energy associated with a scalar field Q 
with self-interaction potential V{Q). If the potential energy density is greater than the kinetic one, then the pressure 
associated to the Q-field is negative 0. Quintessence models that accomodate the present day acceleration tend 
to accelerate eternally, and, as a consequence, the resulting space-times exhibit event horizons ||. A wide class of 
quintessence models, with eternal acceleration, associated with static metrics, have been prezented in ||. 
CN| , Brans-Dicke (BD) theory may explain the present accelerated expansion of the Universe without resorting to a 
cosmological constant or quintessence matter . The conditions under which the dynamics of a self-interacting BD 
field can account for the accelerated expansion have been considered in - Jl^ . Accelerated expanding solutions 
can be obtained with a quadratic self-coupling of the BD field and a negative coupling constant u) A cosmic 

| \ fluid obeying a perfect fluid type equation of state cannot support the acceleration |l3| . The nature of the scalar field 
&Tj. potential compatible with a power law expansion in a self-interacting BD cosmology with a perfect fluid background has 
been analyzed in (lj] . Models with non-minimal coupling of the scalar field have been considered in . Complicated 
cosmological scenarios, with a four-dimensional effective action connected with supergravity and string theory, have 
^ ' been obtained in fL6fl . 

It is the purpose of the present Letter to consider some exact classes of solutions of the field equations in the 
framework of BD theory with a quartic potential, V(</)) ~ </> 4 . By means of some appropriate transformations, the 
field equations can be reduced to a system of two independent Riccati's type differential equations. Three classes of 
exact solutions of the field equations are presented, and their physical properties are investigated in detail. 



II. FIELD EQUATIONS, GEOMETRY AND CONSEQUENCES 



The physical model we are considering is the Brans-Dicke action, along with a self interacting potential V (<fi) , 
coupled to the matter field Lag rangian Lm via the action 



s = 



< f >R -^o t ^ a _ v ^ +Lr 



(1) 
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where uj is the BD coupling parameter. For w = — 1 this action is identical to the low energy string theory action. In 
the present Letter we use units so that &t:G = c = 1. 

For a homogeneous flat space-time, with scale factor a, filled with a perfect fluid, with pressure p m and energy 
density p m , the Einstein-Brans-Dicke gravitational field equations are: 



T 2 _ Pm + P</> n fr , orj2 Pm + P<f> 



3H* = rm ' rw , 2H + 3i? 2 = - em ' , (2) 



M" 1 = ^ + ^ - 3^,P^" 1 = fl + ^)^ 2 -^+^ + 22ty, (3) 



where and are the energy density and pressure associated to the BD scalar field, given by 

where we denoted H = da -1 and ip = cpcp^ 1 . ^ gives the rate of change of the gravitational constant G(t), tp = — §^y- 
The wave equation for the BD field takes the form 



v (2a; + 3) (2a; + 3) V ^ # J 



(4) 



The energy conservation of the matter implies p m + 3 {p. m + p m ) H = 0. 

We assume that the self interaction potential V(<p) has the form V ((f)) — Vo</> 4 , Vq = constant. This form corresponds 
to chaotic inflation model, where the constant Vo is subject to the constraint Vo < 10 -12 , coming from the observational 
limits on the amplitude of fluctuations in the cosmic microwave background |T^ ] . 

With this choice, and with the use of Eqs.(||)-(||), we obtain the following equation describing the dynamics of the 
Universe: 



tp--H = -H 2 - -i/j 2 - 3Htp. (5) 
uj uj 2 



We introduce the following matrices 
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L= [1>--H \,M 1 = (</> H),A=[ 2 J (6) 
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Then Eq.(|) can be written as a matrix equation in the form L = M T AM . 

We define two new variables tp, H , which are obtained by means of a linear transformation, described by the matrix 
K, and applied to the matrix M 

where / + , /_, K + and K- are real numbers. 

In the new variables \tp,H) we have L — M T K T AKM. We choose the elements of the matrix K so that K T AK 



is a diagonal matrix, with the the diagonal elements given by the eigenvalues of A. The eigenvectors and the 
eigenvalues of the matrix A are given by ( f+KZ 1 KZ 1 ) T , ( f-KZ 1 KZ 1 ) T where X T = 12 -^ s , f T = 12 +^ s , 
Kzf = 4^/^ + 1 and s = \/\AA + 24a; + 37a; 2 . s satisfies the condition s > 0,Vo;. 
Therefore the equations satisfied by the new unknown functions ip, H arc 

a_^-A_^ 2 = a+^ + A + H 2 = F(i), (8) 
where we have introduced a solution generating function .F(i) and denoted a^p = ei-fQl 1 (f± — 3a; -1 ), with e± = ±1. 
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III. CLASSES OF EXACT SOLUTIONS OF THE FIELD EQUATIONS 



In order to obtain some classes of general solutions of Eq. (||), we assume first that i[> and H are given by 
■0 = c_< _1 , H = c+t^ 1 , c± = const. These functional forms of the new variables satisfy Eq. if the consistency 
condition a_ — a + , relating the parameters a^, and A^ holds, with = — a^&p + s^zX^c?^. By choosing 
F(t) = a-t~ 2 = a+t~ 2 , we separate Eq. (m into a system of Riccati's differential equations given as 



dip 



dH 



<z_-f - = — ,a+— + X+H 2 = -f 



dt 



t 2 



dt 



t 2 



(9) 



The mathematical form of F(t) has been chosen in order to obtain an exact solution of Eqs. (||). The system (Q) has 
the particular solutions 0o = c_i _1 and Ho = c + t _1 . With the help of the standard transformations -0 = it -1 + ipo, 
H = v^ 1 + Hq the Riccati's differential equations are transformed into two linear Bernoulli's equations. 

Therefore a first class of general solutions of Eqs. (§) is given by 



c_(l+2n_) 



t 2r, 



— +c_t- 1 ,H 



t -2n+ 



d + + c + (l-2n+) t 2 " + 



~ + c+t-\ 



where we have denoted = c^A^a-p 1 . d^ are arbitrary constants of integration. 

With the use of the linear transformation (m) the solution for (■*/>, H) is obtained in the form 



f+KV 1 ^ + f-K- l c + , 
= p + UK+ 1 - 



c_(l+2n_) 



t 2 



rl _1_ "+ f-2w + + l ' 
"+ + c+h-2n + ) t + 



(10) 



(11) 



K^c- + KZ 1 c A 
t 



K 



-l 



2n_ 



t -2n+ 



c_(l+2n_) ' 

The BD scalar field and the scale factor of the Universe are: 

f+ K +' 



+ c + (l-2n+) 



-t 



c- (1 + 2n_) 



c+ (1 - 2n+) 



_^-2n + +l 



(12) 



(13) 



a = a t + - + 



c_ (1 + 2n_) 



2n_+l 



c+ (1 - 2n+) 



2n+ + l 



(14) 



where ao > and (f>o > are constants of integration. 

Another class of exact solutions of the gravitational field equations in the BD theory with matter fluid is obtained 
by assuming for the generating function the form F(t) = (3 = const.. In this case equations (Rh take the form 



a_^ = \^ 2 +P ia+ ^- = -\ + H 2 + r3. 



(15) 



Depending on the sign of the parameters A_, A + and (3, we obtain two distinct classes of solutions. In the first case 
we assume that A_ > 0, — A + > and (3 < 0. With this choice the general solutions of Eqs. (151) are given by 



1> 



XZ 1 coth [AT_ (t - t a )} , H = -J- \/3\ XZ 1 coth [N+ (t - t )] 



(16) 



where N± — y/e^~X±\j3\a± ■ This solution is mathematically consistent for values of -0 and H so that tp > V> < 



'P and U 



ML 



With the use of the transformation (R) we obtain 

- -./+X+y^|A: 1 coth [JV_ (t - t Q )} - f-KZ 1 ^- \/3\ A; 1 coth [N+ (t - to)] 



(17) 
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and 



H = -K^-J\P\Xl l "coth [AG {t-to)} -KZ l \j- |/?| A+ 1 coth [A+ (i-t )] 
On integration, we obtain the BD scalar field and the scale factor: 



f+ K + 1 \/^\ X - 1 l-K 1 ^J-\(3\\ + 1 

b sinh N - [JV_ (t - t )} sinh N + [N+ (t - to)] , 



(18) 



(19) 



a = a sinh 



[AG (t - t )} sinh 



[N+ (t - t Q )} 



(20) 



A third class of solutions is obtained by assuming A_ > 0, — A+ > and > 0. The solution can be formally 
obtained from the previous one by means of the substitution (3 — > —[3, cothix — * i cotx etc. In all these cases the 
energy density and pressure of the baryonic matter and of the scalar field follow from the field equations (||)-(||). 



IV. DISCUSSIONS AND FINAL REMARKS 



The first class of solutions depends on the set of five arbitrary constants c_, d±, Vq and u>. To determine them 
from the actual observational data we take the matter pressure to be zero, p m = 0. For the age of the Universe and 
the Hubble constant we adopt the values t — lAGyr and Ho — Qbkms^ 1 Mpc^ 1 [Q. The density parameters of the 
matter and of the BD scalar fields are defined as f2 m = 3 ^g^ = 0.25 and = jjprj; = 0.75, respectively Q|. For the 

deceleration parameter q — 4i-jf — 1, which is an indicator of the accelerating behavior, we assume an actual value of 
q = —0.5. Therefore, by taking uj as a free parameter, we have four observational constraints to be satisfied by the 
model. Thus the numerical values of the constants c_, d± and Vq can be obtained from fitting the model with the 
observations. Generally, a physical solution for the resulting non-linear system of algebraic equations can be obtained 
only for small negative values of u>. 

The variation of the energy density of the matter and of the BD field is represented, for the first class of solutions, 
in Fig. 1. Generally, the energy density of the scalar field dominates the matter energy density, thus providing the 
dominant contribution to the total energy density. The deceleration parameter is represented, for different values of 
lu < 0, in Fig. 2. q is negative, with values in the range — 1 < q < 0, indicating an accelerating evolution of the 
Universe. 
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FIG. 1. Time evolution of the energy density of the matter p m (solid curve) and of the scalar field p$ (dashed curve) (in 
units of 10~ i7 GeV 4 ) for the first class of solutions for u> = -2 ,d_ = 0.32, d+ = 0.53, c_ = -1.95, Vo = 0.0064 and c+ = -0.54. 
The time is expressed in units of 10 Gyr. 
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FIG. 2. Dynamics of the deceleration parameter for the first class of solutions for lo = —1.8 (solid curve), for lo = —2 (dotted 
curve), and for u> — —2.2 (dashed curve). We have used the values d_ = 0.32, d + = 0.53, c_ = —1.95, Vq = 0.0064 and 
c+ = —0.54. The time is expressed in units of 10 Gyr. 

The second class of solutions depends (by chosing to — 0) on three parameters (3, Vq and lo. Their numerical values 
can be obtained so that the solution fits the actual observational values of H, £l m and 0^,. The resulting algebraic 
system of equations has also a physical solution only for small negative values of lo. 

The time variations of the energy density of the matter and of the BD field and of q for the second class of solutions 
are represented in Figs. 3 and 4. The energy density of the scalar field dominates the matter energy density, and the 
evolution is accelerating. 
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FIG. 3. Time evolution of the energy density of the matter p m (solid curve) and of the scalar field (dashed curve) (in 
units of 10 _47 GeV 4 ) for the second class of solutions, for lo = —2, [3 = 0.001 and Vb = 0.0073. The time is expressed in units 
of 10 Gyr. 
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FIG. 4. Dynamics of the deceleration parameter in the second class of solutions for u = — 1.8 (solid curve), lo = — 2 (dotted 
curve) and lo — —2.2 (dashed curve). In all cases we have used the values /3 = 0.001 and Vb = 0.0073. The time is expressed 
in units of 10 Gyr. 



Therefore, these two classes of solutions are compatible with the observed cosmological data only for small negative 
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values of the coupling parameter of the BD field. An other important observational bound requires fi^, < 0.044 at 
nucleosynthesis [|l 8f| . This condition can also be satisfied by appropriately chosing the matter equation of state in the 
very early stages of evolution of the Universe. Generally, in the present models we have > p m , and the evolution 
is accelerating. An accelerating expansion during the whole evolution of the Universe is not favored by the curent 
models of structure formation. However, as shown in |19|, the scalar field, which is non-minimally coupled to gravity, 
may undergo clustering processes, eventually forming density perturbations, which can be investigated only within 
a non-linear approach, and leading to the formation of some cosmic structures. On the other hand, to satisfy the 
actual observational constraints, a much larger value than that predicted by the inflationary scenario is needed for 
the constant Vq. 
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